Contributions of hyperon-hyperon scattering to subthreshold cascade production in 

heavy ion collisions 
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Using a gauged flavor SU(3)-invariant hadronic Lagrangian, we calculate the cross sections for 
the strangeness-exchange reactions YY -f-^- A''H {Y — A, E) in the Born approximation. These cross 
sections are then used in the Relativistic Vlasov-Uehling-Uhlenbeck (RVUU) transport model to 
study H production in Ar-fKCl collisions at incident energy of 1.76A GeV and impact parameter 
b = 3.5 fm. We find that including the contributions of hyperon-hyperon scattering channels strongly 
enhances the yield of H, leading to the abundance ratio E~/(A -I- E") = 3.38 x 10"'', which is 
essentially consistent with the recently measured value of (5.6 ± 1.2j;j'7) x 10~^ by the HADES 
collaboration at GSI. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 



The study of particle production in heavy ion colhsions 
at energies below their thresholds in nucleon-nucleon 
collisions was a topic of extensive studies during the 
1990s [ll~lH|- The main motivation for such study is that 
it offers the possibility of extracting information on the 
nuclear equation of state (EOS) at densities above that of 
normal nuclear matter. In particular, the yield of strange 
hadrons, such as the kaon, has been shown to be sensi- 
tive to the stiffness of the nuclear equation of state up to 
three times normal nuclear matter density, with a softer 
EOS giving a larger yield than a stiff EOS. Indeed, exper- 
imental results obtained by the KaoS Collaboration Q 
at the Society for Heavy Ion Research (GSI) in Germany 
on the yield of kaons in heavy ion collisions at subthresh- 
old energies have led to the conclusion that the nuclear 
equation of state at high densities is soft, consistent with 
an incompressibility of about 200 MeV extracted from 
the collective flow studies by the Plastic Ball and 
EOS [1] Collaborations from Lawrence Berkeley Labo- 
ratory (LBL) and the E877 [I and E895 ^ Collabora- 
tions at the Alternating Gradient Synchrotron (AGS) of 
Brookhaven National Laboratory (BNL). More recently, 
the doubly strange baryons S from Ar-fKCl collisions 
at 1.76 A GcV, which is below the threshold energy of 
3.74 GeV in a nucleon-nucleon collision, was measured 
by the HADES Collaboration at GSI The mea- 

sured abundance ratio including the statistical and sys- 
tematic errors is S-/(A -\- S°) = (5.6 ± 1.2;^;^) x lO'^. 
This value is about 10-20 times larger than those given 
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by the statistical model and the relativistic trans- 
port model [l3|- Because of the very low collision en- 
ergy, secondary reactions other than the direct reaction 
NN — )■ N'E.KK are expected to contribute significan- 
tantly to S production in these collisions. In Ref. [l3j . 
the strangeness-exchange reaction KY — ?> ttS {Y = A, E) 
between antikaon and hyperon was introduced in the 
Vlasov-Uheling-Uhlenbeck (RVUU) transport model [l3l 
to study 5 production in heavy ion collisions. The cross 
sections used in Ref. [l^l were taken from the coupled- 
channel calculation of Ref. [l5| based on a gauged fla- 
vor SU(3)-invariant hadronic Lagrangian. Since there 
are more hyperons than anitkaons in heavy ion colli- 
sions at this energy, the strangeness-exchange reaction 
YY — >■ iVS between two hyperons is expected to be 
important for S production in these collisions. In the 
present study, we use the same hadronic Lagarangian 



as in Ref. 
reaction YY 
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to evaluate the cross sections for the 
iVS. For an exploratory study, these 



lwchen@sjtu.cdu.cn 
ko@comp.tamu.edu 
,suhoung@yonsei.ac 



cross sections are calculated in the Born approximation 
with the cutoff parameter in the form factors at interac- 
tion vertices fitted to the cross sections for the reactions 
KY — ^ ttS obtained in Ref. [l^. For completeness, we 
also include the reaction KN — >■ K'E. with its cross sec- 
tion taken from empirically available values. Our results 
show that the inclusion of the reaction YY — > N'E. sig- 
nificantly enhances the yield of S in heavy ion collisions 
at subthreshold energies, resulting in the abundance ra- 
tio S-/(A -I- E") = 3.38 X 10-3 in Ar-fKCl collisions at 
1.76 A GeV and impact parameter & = 3.5 fm, which is 
essentially consistent with the recently measured experi- 
mental value. We find, however, that the contribution of 
the reaction KN — >■ K'E. to the S yield is negligible. 

The paper is organized as follows. In Sec. II, we 
describe the gauged flavor SU(3)-invariant hadronic La- 
grangian [l5| . calculate the amplitudes for the reaction 
YY N'E in the Born approximation, and parametrize 
the resulting cross sections. In Sec. HI, we introduce 
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the paramctrization of the empirical cross section for the 
reaction KN KE. as a function of the center of mass 
energy. Wc then briefly review in Sec. IV the RVUU 
transport model for high energy heavy ion collisions. Nu- 
merical results on the time evolution of the S abundance 
in Ar+KCl collisions at 1.76 A GeV and impact parame- 
ter b — 3.5 fm arc presented in Sec. V. Finally, wc present 
some discussions in Sect. VI and a summary in Sec. VII. 

II. THE HADRONIC MODEL 



they are included by replacing the partial derivative 9^ in 
Eq.dl]) with the covariant derivative = 9^ — [V^,], 
where g is another coupling constant. 

We further include the tensor interactions between 
baryons and vector mesons via the interaction La- 
grangian 

C} = I^Tr [(2a - 1) Ba^^Bd^y, + Bo^'' B\ , 

(5) 

with being the tensor coupling constant. 



Possible reactions for S production from hyperon- 
hyperon collisions arc AA — > A^S, AE — > A^S, and 
EE — )■ A^S. Cross sections for these reactions can be eval- 
uated using the same Lagrangian introduced in Ref 
for studying 5 production from the reactions Kh. — >■ ttS 
and A'S — > ttS. This Lagrangian is based on the gauged 
SU(3) flavor symmetry but with empirical masses. The 
coupling constants are taken, if possible, from empirical 
information. Otherwise, the SU(3) relations are used to 
relate unknown coupling constants to known ones. Also, 
form factors are introduced at interaction vertices to take 
into account the finite size of hadrons. 



A. The Lagrangian 

As in Ref. [l5|, we use the following flavor SU(3)- 
invariant hadronic Lagrangian for pseudoscalar mesons 
and baryons 

L = iTT{BpB) +T:r[{d^P^d''P)] 

+ g' {Tr [{2a - 1) Bj^^Bd^P + 57^7" (d^P) B] } , 

(1) 

where B and P denote, respectively, the baryon and pseu- 
doscalar meson octets 
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with g' being a coupling constant and a being a param- 
eter. 

For the interactions of baryons and pseudoscalar 
mesons with the vector meson octet Vu, 



V 




K* 



(4) 



B. Born approximation to the reactions AA N3, 

AE NE, and EE ^ NE 



-N A. E 




FIG. 1: Born diagrams for the reactions A A — >■ AH, AE — >■ 
NE, and EE ^ A^H. 

In the Born approximation, the reactions AA — >■ N'E., 
AE AS, and EE — ?> A^S are described by the tree- 
level i-channel and u-channel diagrams shown in Fig. 1. 
To evaluate their amplitudes requires the following inter- 
action Lagrangian densities that are deduced from the 
hadronic Lagrangian in the previous subsection, i.e., 

Ckna = l^NJ^^^^A^f,K + U.c., 
rriK 

Cknj: = :^i^A757^(f.E)a^i^ + H.c., 
rriK 

Ckae = ^S7 VAa^A'^ + H.c, 
niK 

Ckse = ^S7V(^-S)9^i^^ + H.c, 
rriK 

Ck'NA = gx'NAN (rK + ^^'^^ af'^d^K: ) A 



rriN + rriA 



+ H.c, 



X E • I Y'rK* + '"^•^^ a'^'^TdJC ) A^ 



iTiN + m-s 



H.c, 



+ H.c, 

C-K'T.E. ~ gK'Y:3 



X E • j^tk;^ + 
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H.c. 



(6) 



In the above, r are Pauli matrices; tt, p, and E denote 
the pion, rho meson, and sigma hyperon isospin triplets, 
respectively; K = {K+,K°Y [K* = {K*+,K*°)^) and 
K" = {K^,~K-y [K*" = {K*°,-K*-)^) denote the 
pseudoscalar (vector) kaon and antikaon isospin doublets, 
respectively; and S = is the cascade hyperon 

isospin doublet. The coupling constants in above interac- 
tion Lagrangian densities are relate to those in Sec. Ill Al 
by 



Jkna 


2a- 3 , 


Ik NT, 


2a - 1 




2^3 ^' 


rriK 


2 


fxAE 


3-4a , 


/ifSH 


1 , 


rriK 


2x/3 ■'^ ' 


ITlK 








x/3 




9k'NK 


= -gK'AE = 


^' 








9 






= .gK*sH = - 


^4' 





KK'AN 



KK'AS 



9k'An 

9k* AN ' 

9k'As 



-, KK'NS 



9k*na 
ruN + TUA 

9k'Ae 
rriA + ms 



9K-A~ 

2a -3 5* 



I^K'Y.E. 



__ 9k'NT. 

gK'NT. ' 
9k'SS 

9k'te ' 



9 k* NT 



2\/3 2m' niN + 
3 - 4a 5* 9k'TS 



2^3 2m' ms 



2a - 1 g* 
2 2^' 

4m 



(7) 



X [E (p,) ^^rYi (pi)] [N (p4) l^'l'-Y^ (P2)] 

+9k'YxB.9k'NY2 

t ~ m 
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with u = (pi —pi)^. The form factor F introduced at 
the interaction vertex because of the hardron structure 
is taken to have the monopole form, 
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and depends on the three momentum transfer q and the 
parameter A. The isospin factors rju, Tjtu = i^ut, and 
rjuu in Eq. (|9]), which are obtained from summing the 
isospins of initial and final particle, are 18, 10, and 18 
for the reaction EN; 6, 2, and 6 for the reaction 

AS EN, and 2, 2, and 2 for the reaction AA EN. 



C. Cross sections for the reactions AA — > A'^S, 
AE AH, and EE A^H 



The cross sections for these reactions are then given by 

1 



CTYY^Nsis) 



64:Trspf 



dt\Mr, 



(8) 



where s ~ {pi +^2)^ and t — {pi — pa)^ arc the usual 
squared center of mass energy of colliding hyperons and 
the squared four momentum transfer in the reaction; and 
Pi is the momentum of initial hyperons in their center of 

mass frame. The spin-isospin averaged amplitude \A4\^ 
in the above equation is given by 



(2si + 1) (2S2 + 1) (2/1 + 1) (2/2 + 1) 



(9) 



TABLE I: Coupling constants used in the present study. 
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For numerical calculations of the cross sections, we use 
the coupling constants shown in Table I. These values are 
obtained from g' = UA GeV"\ g = 13.0, and gy2m = 
19.8/mN that are determined from the empirical values 
Unn = 1-00, gpNN = 3.25, g*^^ = 19.8 [ll], and a = 

0. 64 using relations based on the SU{3) symmetry, 

1. e., 



where M* ^ , , and Af " ^ , , arc the spin-dependent 
amplitudes for the two Born diagrams shown in Fig. [1] 
and are given by 
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For the cutoff parameter A in the form factor, its value 
is taken to be A = 0.7 GeV in order to reproduce, 
as shown in Fig. [51 the cross sections for the reactions 
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FIG. 2: (Color online) Isospin-averaged cross sections for (a) 
KA — >■ ttH and (b) KT, — >■ ttH. Solid lines are from the Born 
approximation with the cutoff parameter A = 0.7GeV in the 
form factor, and dashed lines are those based on the coupled- 
channel calculation [l^ . [isj . 



KA — >■ 7r5 and KI] -> ttS that are obtained from the 
coupled-channel calculation based on the same hadronic 
Lagrangian [l5| . 

In Fig. 121 we show by solid lines the isospin-averaged 
cross sections for the reactions A A NE (panel (a)), 
AE NE (panel (b)), and ES NE (panel (c)) as 
functions of the center-of-mass energy yG, obtained with 
A = 0.7 GcV. These cross sections can be parametrized 
as 
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where pA and pN are initial A and final nucleon momenta 
in the center-of-mass frame. We note that the magnitude 
of our cross section for the reaction EN — > AA is similar 
to that of Ref. [l^ obtained from the SUg quark model 
formulated in the resonance group method but is smaller 
than that extracted from the {K^ , K'^)E~ reactions in a 
nucleus [l^. For comparisons, we also show in Fig.|3]the 
cross sections for the reaction YY — )• NE for the cutoff 
parameters A ~ 0.5 GeV (dashed lines) and A = 1 GeV 
(dotted lines). As expected, the cross sections are larger 
for a larger A. 

The cross sections for the inverse reactions (Jne^aa, 
CTArH^-Ai;, and ctateh-ss are related to above cross sections 
by the detailed balance relations: 
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FIG. 3: (Color onhne) Cross sections for (a) AA — > A^H, (b) 
AE ^ A^H, (c) EE ^ A'H, (d) A^H ^ AA, (e) A^H ^ AE, 
and (f) NE, — > EE as functions of the center-of-mass energy 
^/s from the Born approximation with cutoff parameters A = 
0.5 GeV (dashed lines), A = 0.7 GeV (solid lines), and A = 
1 GeV (dotted lines). 



III. CROSS SECTIONS FOR THE REACTION 

KN KE 



For completeness, we also include in the present study 
the reaction KN — KE. Both the differential and total 
cross sections for this reaction were measured in 1960s 
and 70s [20l - |29j , and they are shown in Fig. by solid 
squares for K~+p K^E" (panel (a)), K~+p — )• A'°S° 
(panel (b)), and K^ + n ^ K'^E'^ (panel (c)). Recently, 
a phenomenological model was introduced in Ref. [soj to 
describe these reactions, and the results are shown by 
dashed lines in Fig. |4l In the present study, we use the 
following parametrization for these cross sections: 
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FIG. 4: (Color online) Cross sections for the reaction K~ + 
p — >■ K^'ST (upper left window), +p — > K'^'EP (upper right 
window), and K~ + n — >■ K^'SP (lower left window). Solid 
squares are experimental data, dashed lines are theoretical 
results of Ref. [s^l, and solid lines are our parametrization. 
Lower right window shows the isospin averaged cross sections 
for the reaction KN — >■ KE. (solid line) and its inverse reaction 
KE. — KN (dashed line) based on our parametrization. 



In terms of these cross sections, the isospin averaged cross 
section for the reaction KN — /-CS can be expressed as 

+ crK-n^K"3-)- (17) 

The detailed balance relation then allows us to express 
the cross section for the inverse reaction KE — > KN as 



Pn 

PS 



(18) 



where pN and ps. arc the 3-momenta of nuclcon and :^ in 
the ccnter-of-mass frame. 



IV. THE RELATIVISTIC 
VLASOV-UHLING-ULENBECK TRANSPORT 
MODEL 

To study S production in heavy ion collisions at sub- 
threshold energies, we generalize the RVUU transport 
model [3l to include the reactions YY — >■ NS and 
KN — >■ if S and their inverse reactions besides the re- 
action KY — > ttS and its inverse reaction that were al- 
ready included in Ref. • In addition to these reactions 
and other reactions involving nuclcons. Delta resonances, 
hyperons, pions, kaons, and antikaons, the VUU model 
also includes the mean-field effect on the propagation of 
baryons, kaons, and antikaons. For nucleons and Delta 
resonances, their mean- field potentials are taken from the 
relativistic mean-field model via the scalar and vector 



potentials, so their motions are given by the following 
equations of motion: 



P_ 

E* 



Wo) 



(19) 



where m* = 7ti — $, p* = p - W, E* = y'p*^ + m*^ 
with $ and W = (VFo,W) being the scalar and vector 
mean fields, respectively. These mean fields are calcu- 
lated from the effective chiral Lagrangian of Ref. [3l| with 
parameters determined from fitting the nuclear matter 
incompressibility Kq = 194MeV and the nucleon effec- 
tive mass m* /m = 0.6 at normal nuclear matter density 
Pq = 0.15 fm'^. For A and E hyperons, their mean-field 
potentials are taken to be 2/3 of the nucleon mean-field 
potential according to their light quark content. Simi- 
larly, the mean- field potential for S is 1/3 of that of the 
nucleon . 

For kaons and antikaons, their mean-field potentials 
are derived, on the other hand, from the dispersion rela- 
tion obtained in the chiral Lagrangian [32| 
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-,1/2 
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(20) 

where ps = (NN) is the scalar density, / = 103 MeV is 
the pion decay constant, and the ± is taken as for 
kaons and "-" for antikaons. The KN and KN sigma 
term S jvx in the above equation can in principle be cal- 
culated from the SU{3)l x SU{3) r chiral Lagrangian but 
are taken to have the values Y^nkI P = 0.22 GeV^fm'^ 
and 'S'nkIP = 0-35 GeV^fm^ as in Ref.fll from fitting 
the kaon and antikaon yields in heavy ion collisions. 

Besides affecting the propagation of particles, the 
mean-field potential also has effect on the threshold en- 
ergy for particle production as a result of the potential 
difference between the initial and final states of a reac- 
tion. For example, this effect is important for under- 
standing the enhanced production of antikaon through 
the reactions BB ^ BBKK, nB o KKB, and nY o 
KN in heavy ion collisions at subthresold energies. As a 
result, the contribution of the reaction KY — >■ ttS to S 
production in heavy ion collisions at subthreshold ener- 
gies was found in Ref. [ll] to be further enhanced. We 
note that in the RVUU model, kaons, antikaons, hyperons 
(lambdas and sigmas), and cascade particles are treated 
perturbatively by neglecting the effect of their produc- 
tion and annihilation on the collision dynamics, which 
is dominated by the more abundant nucleons. Delta res- 
onances, and pions. In this approach, kaons, antikaons, 
and hyperson are produced from nucleon (Delta)-nucleon 
(Delta) and pion-nucleon (Delta) collisions whenever it 
is energetically allowed, and they are given probabili- 
ties that are determined by the ratios of their respective 
production cross sections to the total cross sections of 
the colliding particles. For S production from antikaon 
collisions with nucleons or hyperons and from hypcron- 
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hyperon collisions, it is similarly treated but the proba- 
bility of the produced S is reduced by the probabilities 
of colliding particles. The annihilation of these rare par- 
ticles is treated in a similar way and leads to reductions 
of their probabilities. The present approach thus takes 
into account the small probability associated with the 
production of two rare particles in a subthrehold heavy 
ion collision that are involved in the production of a S. 



RESULTS 
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FIG. 5: (Color online) Time evolutions of (a) central baryon 
density (right scale) and the abundances (left scales) of tt, 
A; (b) K, A, E, and K\ and (c) H produced from different 
reactions. 



(right scale). It is seen that the colliding system reaches 
its highest density of about 1.87po at about 7 fm/c when 
most particles are produced. The tt abundance saturates 
at 10.3. Assuming isospin symmetry, the 7r~ number is 
then 3.43 which is very close to the measured number 
of 3.9 ± 0.1 ± 0.1 by the HADES Collaboration [H, 113. 
The time evolution for the abundances of K, A, and S 
arc shown in Fig.[5ljb), and they saturate at the values of 
5.32 X 10-2, 1.15 X 10-3, 2.60 x IQ-^, and 2.60 x 10"^ re- 
spectively. Assuming isospin symmetry gives 2.61 x 10-^ 
for the number, 5.75 x IQ-'' for the K~ number, 
and 3.47 x 10-^ for the A -I- E° number. These num- 
bers are again close to corresponding measured numbers 
of (2.8±0.2±0.1±0.1)xl0-2, (7.1±1.5±0.3±0.1)xl0-^ 
and (4.09 ± 0.1 ± 0.17) x 10-^ by the HADES Collabo- 
ration [ssj . For the time evolution of the S abundance, 
it is shown by the solid curve in Figl5l^c) and is seen to 
saturate at the value 2.34 x 10^^. Taking S- as half of S 
by assuming isospin symmetry, we obtained a 2- num- 
ber of 1.17 X IQ-^ which is about half of the measured 
number of (2.3 ± 0.9) x 10-'' by the HADES Collabora- 
tion [3^. Our results thus lead to an abundance ratio 
S-/(A -I- E°) = 3.38 X 10-3, which is essentially consis- 
tent with the measured value of (5.6 ± 1.21^ 7) x 10-^ by 
the HADES collaboration. 

The contributions to S production from different reac- 
tion channels are also shown in FiglS^c). Dotted, dashed- 
dotted, and dash lines denote, respectively, the abun- 
dance of the S particles from the reactions YY — >■ A^S, 
KY — >• ttS, and KN — )- ifS. Compared to the total S 
abundance, shown by the solid line in FiglSl the contri- 
butions are 97.5%, 2.40%, and 0.1% from the reactions 
YY -)■ A^S, KY ttS, and KN -> ATS, respectively. So 
the YY — )■ N'E. channel dominates S production in heavy 
ion collisions at subthreshold energies. This can be ex- 
plained by the fact that the cross section for YY -> A^S 
is almost 3-4 times the cross section for KY — >■ ttS, and 
almost hundred times the cross section for KN — >■ ATS. 
Also, the hyperon abundance in the system is almost 
20 times the anti-kaon abundance. We note that the 
relative contributions to the S yield from the reactions 
AA NE, AE NE and EE NE arc about 1, 4 and 
1. 



In this Section, we show the results for ^^Ar -I- KCl 
collisions at incident energy 1.76 AGeV, taking as an av- 
erage of ^"Ar + K^^ collisions and "'"Ar + Cl"^^ collisions, 
and compare them with the data from the HADES Col- 
laboration at SIS. The HADES trigger (LVLl) selects 
approximatel y th e most central 35% of the total reaction 
cross section [ll| . According to GEANT simulations [s^ 
with the UrQMD [13, [13 transport approach as event 
generator, the average value and width of the correspond- 
ing impact parameter distribution amount to 3.5 and 1.5 
fm, respectively. For simplicity, we take b — 3.5 fm in the 
present study. Fig.[SJa) shows the time evolution of tt and 
A abundances (left scale) and the central baryon density 



VI. DISCUSSIONS 

Our results are obtained without the consideration of 
the isospin asymmetry effect due to different proton and 
neutron numbers in the colliding nuclei, which is ex- 
pected to increase the final abundance ratio S-/(A+E'^). 
If we assume that the abundance of E has reached chem- 
ical equilibrium in heavy ion collisions, which is certainly 
questionable in view of the failure of the statistical model 
in describing the experimental data, this enhancement 
can be estimated using E~ /E° = c^i^^/'^ — /TP = 
E°/E+ = N/Z, where /xc is the charge chemical poten- 
tial and T is the temperature of the system. With the 
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value N/Z ~ 1.14 for Ar'"' + K^^ or Ar"" + Cf^, we have 
S" = 0.533 S and TP = 0.3314 E, leading to the ratio 
I{K + YP) = 3.60 X 10-3 that is 6.5% larger than that 
for an isospin symmetric system. 

Also, the nuclear EOS used in the transport model can 
affect the final S abundance in heavy ion collisions. The 
results presented in the previous Section arc based on a 
soft EOS. Using a stiff EOS, we find that the A, S, and S 
abundances are reduced to 1.74 x 10"^, 1.77 x 10"^, and 
1.46 x lO"'^, respectively. The reason for this reduction 
in the hyperon abundances is that the energy density 
of the colliding system increases faster for a stiff EOS, 
thus making its expansion faster and reaction time short. 
However, the abundance ratio 5"/(A+S'') = 3.13 x 10^"^ 
for the stiff EOS is essentially the same as that for a soft 
EOS. 

Furthermore, the results presented here are for the im- 
pact parameter 6 = 3.5 fm. A more realistic comparison 
with experimental data should include a distribution of 
impact parameters. We have checked that using differ- 
ent impact parameters, the ratio S-/(A -I- YP) remains, 
however, essentially unchanged, since both hyperons and 
cascade abundances change by almost the same factor. 

Finally, because of the very large S production cross 
sections and the small size of the colliding system, the 
geometrical treatment of S production from hyperon- 
hyperon scattering in terms of their scattering cross sec- 
tion as used in the RVUU transport model may become 
inaccurate. This can be seen from the dependence of fi- 
nal S abundance on the value of the cutoff parameter A 
in the form factor used in evaluating the cross sections 
of the reactions YY — > N'E.. As shown in Fig. |3l these 
cross sections increase with increasing value of A. Re- 
sults from our transport model study show, on the other 
hand, that the S abundance increases with decreasing 
value of A. However, our conclusion in the present work 
is expected to remain unchanged since the S abundance 
changes only by about 30% when the S production cross 
sections change by more than a factor of 4. We note that 
a more accurate treatment of particle scattering may be 
achieved by using the stochastic method of Ref. [40| based 
on the transition probability, and we hope to purse such 



an improved study in the future. 



VII. SUMMARY 

We have calculated the cross sections for the reac- 
tion YY N^{Y = A, E) based on a gauged SU(3)- 
invariant hadronic Lagrangian in the Born approximation 
and found that these cross sections are almost 4 times the 
cross sections for the reaction KY — )- ttS that was con- 
sidered in previous studies. We then used these cross sec- 
tions to study S production in ^°Ar-|-KCl collisions at the 
subthreshold energy of 1.76 AGeV within the frame work 
of a relativistic transport model that includes explicitly 
the nucleon. A, pion, and perturbatively the kaon, an- 
tikaon, hyperons and S. We found that the reaction 
YY -> N'E. would enhance the S abundance by a factor of 
about 16 compared to that from the reaction KY — >■ ttS, 
resulting in abundance ratio S^/(A + TP) = 3.38 x lO-'^ 
that is essentially consistent with that measured by the 
HADES Collaboration at GSI. Our study has thus helped 
in resolving one of the puzzles in particle production from 
heavy ion collisions at subthrchold energies. 
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